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Abstract 

A simple method to explore the interaction of antihyperons in nuclei by exclusive hyperon-antihyperon pair production 

plose to threshold in antiproton nucleus interactions is proposed. Due to energy and momentum conservation event- 

OO by-event transverse momentum correlations of the produced hyperons and antihyperons contain information on the 

2^ [difference between their potentials. A schematic Monte Carlo simulation is used to illustrate the sensitivities of the 

Z^ proposed method for the reaction 1.66 GeV/c p^^C -^ AA. For produced D-meson pairs at 6.7 GeV/c the sensitivity of 

'the transverse momenta correlation will probably be too small to deduce differences between the potentials for D^ and 

TT^'D" mesons. However, for SS pairs produced at 2.9 GeV/c the asymmetry is sufficiently sensitive to predicted differences 

ry^ between the 5 and S potentials. 
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Based on G-parity transformation 1] Diirr and Teller 
predicted within an early form of a relativistic field the- 
pry a strongly attractive potential for antiprotons in nu- 
clei [^ y] . It is however obvious that G-parity transforma- 
I Ition can provide a link between the NN and NN interac- 
tion at most for distances where meson exchange is a valid 
[concept y, y . For distances lower than about 1 fm, quark 
degrees of freedom may play a decisive role. The study of 
,the potential of antibaryons in nuclei may therefore help to 
elucidate the role of the quark-gluon structure of baryons 
,for the short-range baryon-baryon force. 

Earlystudies of antiproton- nucleus scattering cross sec- 
jtions [a, lZI showed however disagreement with such a strong 
attractive potential. Later, X-ray transitions in antipro- 
'tonic atoms [8|, |9|, [lOl, [ill (an overview on subsequent ex- 
perimental studies can be found in Ref. [ij]) gave also 
,hints for an attractive potential albeit with large uncer- 
tainties [iSl, [141 • Other analyses favor shallow real and 
[deep imaginary potentials (for example [15|). More com- 
prehensive studies [16| of antiprotonic X-rays as well as 
recent analyses of the production of antiprotons in reac- 
tions with heavy ions resulted in real attractive potentials 
in the range of about -100 to -150 MeV [HHEl. 

Concerning baryons beyond SU(2), only for A hyper- 
ons reliable information on their nuclear potential is avail- 
able from hypcrnuclci studies. No experimental informa- 
tion on the nuclear potential of antihyperons exists so 
far. Mishustin and co-workers recently suggested to study 
deeply bound antibaryonic nuclei via various characteris- 
tic signals in their decay process [20|, l21| , like the produc- 
tion of multi-quark-antiquark clusters, multifragmentation 



events with strong radial flow or sharp lines in meson spec- 
tra due to transitions from the Fermi to the Dirac sea. 
From the experimental point of view it is however not 
obvious whether these proposed observables will provide 
unique and quantitative signals of deeply bound antibary- 
onic systems. 

In this letter we show that quantitative information 
on the antihyperon potentials relative to that of the corre- 
sponding hyperon may be obtained via exclusive antihyper- 
on-hyperon pairs production close to threshold after an 
antiproton-proton annihilation within a complex nucleus 
(Fig. [T|). Once these hyperons leave the nucleus and are 
detected, their asymptotic momentum distributions will 
reflect the depth of the respective potentials. A deep po- 
tential for one species could result in a momentum dis- 
tribution of antihyperons which differs from that of the 
coincident hyperon. This situation is in line with the case 
of antiprotons produced in heavy ion collisions close to 
threshold [17|, [jjl, Il9|. The advantage here is, that we 
are dealing with a quasi stationary system having a rea- 
sonably well defined geometry and that the kinematics is 
determined essentially by energy and momentum conser- 
vation of a (nearly) two-body reaction. However, since 
in the pp center-of-mass the distribution of the produced 
baryon-antibaryon pair will usually not be isotropic, the 
analysis can rely only on the transverse momenta of the 
outgoing baryons: because the initial average transverse 
momentum is equal to zero and neglecting for the moment 
rescattering and the Fermi motion of the struck proton 
(both effects will be discussed below) the magnitude of 
the transverse momenta of the produced baryons and an- 
tibaryons will be equal unless there is a difference in the 
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Table 1: Scalar and vector potentials, S and V, used in the model calculations. The first three columns give the default values for p p 
and A hyperons used in the parameter scans which are presented in Fig. |2] and Fig. [3] The values in the last 8 columns were adopted from 
Ref. [38l . 1 39( 1 and are used in the calculations shown in Fig. 141 
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Figure 1: Scheme of the reaction proposed to explore the nuclear 
potentials of the baryon and the antibaryon (top). The lower part 
shows the probability density distribution of the impact parameter 
leading to the observation of a AA pair in 1.66 GeV/c p^^C— > AA 
reactions. In these simulations default parameters as discussed in 
the text have been used. For orientation the dashed line gives the 
assumed radial density of the ^^C nucleus (right scale). 



effective potentials. 

In the following we explore the influence of the po- 
tentials on the transverse momentum distributions of the 
coincident hyperons and antihyperons as well as on their 
event-by-event correlations by means of a schematic Monte 
Carlo simulation. Albeit crude, this classical approach al- 
lows to explore the role of different features of the reac- 
tion in a transparent way. As an example we consider 
the p^^C — > AA reaction at 1.66 GeV/c, where existing 
data 23, |23| demonstrate the feasibility of such measure- 
ments. Since the method relies essentially on momentum 
and energy conservation, an extension to other hadron- 
antihadron pairs produced exclusively in antiproton-nucleus 
collisions close to their respective thresholds is straight for- 
ward. In future, such reactions can be studied at the inter- 
national Facility for Antiproton and Ion Research FAIR 
24| with e.g. the planned PANDA experiment |25|. 

The absorption of the antiprotons entering the target 



nucleus determines the points of annihilation inside the 
nucleus and the paths which the eventually produced hy- 
perons and antihyperons have to pass inside the nucleus 
prior to emission. For the proton density we adopted a 
Fermi-type distribution 

Po 



p{r) 



(1) 



l + exp{{r-rQ-A^/^)/do)' 

where r denotes the radial distance from the center of the 
target nucleus and A its mass number. For the radius pa- 
rameter tq and the surface diffuseness do default values of 
1.07 fm and 0.54 fm were used [26]. Since in the following 
we focus on light nuclei, equal neutron and proton density 
distributions were assumed. 

The initial pp annihilation is controlled by an pN an- 
nihilation cross section of 50 mb [27|. The lower part of 
Fig. [1] shows the probability density distribution of the 
impact parameter leading to the emission of a AA pair 
in 1.66 GeV/c p^^C -^ AA reactions. For orientation the 
dashed line gives the assumed radial density profile of the 
-•^^C nucleus. Because of the strong absorption of the an- 
tihyperons, the emitted hyperon-antihyperon pairs are - 
unlike in inclusive reactions 28, 29l - created close to the 



corona of the target nucleus at an average impact param- 
eter of 3.1 fm and a typical density of 20 to 25% of the 
central nuclear density. 

For both, the emitted A's and the A's the inverse of the 
average integrated path weighted with the local density 
along the path {j pds)^^ varies for the parameter range 
discussed in this paper in the range from about 800 to 1600 
mb. For our default parameters this value is about lOOOmb 
and thus significantly larger than the typical elastic cross 
sections in the relevant momentum range of <200mb. As a 
consequence re-scattering effects are expected to be small 
and have been neglected in our model. Experimentally, 
re-scattering effects with momentum transfers beyond the 
typical Fermi momentum can possibly be reduced by con- 
straining the azimuthal angle between the hadron and an- 
tihadron momentum. 

In reactions close to threshold the Fermi motion of the 
protons inside the nuclear target contributes significantly 
to the final momenta. Hence the initial proton momentum 
was sampled from a distribution 



dP(p,0,^)oc(l+e(^-^-)/'^-^)- 



V^ sm{e)dp de dcf) (2) 



where a default Fermi energy Ep corresponding to a mo- 
mentum oi pp= 220 MeV/c was used [30[. Quasifree me- 
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Figure 2: Transverse momentum of A hyperons (open symbols) and A hyperons as a function of the total A potential for different parameter 
sets. The dashed and solid arrows mark the result for A and A if all individual A and A potentials are set to 0. 



son scattering experiments suggest [31|, |32|, |33| that owing 
to the fact that the pp annihilations happens in the periph- 
ery of the target nucleus at subsaturation density, signifi- 
cantly lower Fermi momenta may be expected. We there- 
fore varied pp in the range of 180 to 260 MeV/c. For the 
diffuseness parameter a value of T= 1 MeV was used and 
the maximum possible energy was determined by Ep+Es, 
where Eb~8 MeV denotes the typical nucleon binding en- 
ergy. 

Antilambdas produced in pp annihilations are emitted 
preferentially in the direction of the incident antiproton 
^22]. For the 1.66 GeV/c p^^C -^ A A reaction [22, 34, 35|] 
the probability distribution of the center of mass angle 6cm 
of the outgoing antihyperon can be described by: 

exp((cos6'„„ - l)/&o) 



dP{9cm) oc 



flo 



-d{cose„n) (3) 



1 + ao 

with default values for the constant term ao and the width 
6o of 0.1 and 0.3, respectively. 

Lacking any detailed experimental information it is 
plausible to assume that the annihilation cross sections 
for antihyperons show a similar momentum dependence as 
the pp system [36|. We therefore parameterized the ab- 
sorption cross section of the A as 



AAf 



100m6/(p^+l) 



(4) 



with the A momentum given in GeV/c. For the A hy- 
perons a momentum independent inelastic cross section of 



20 mb was adopted. 

The energy and the momentum of the baryons propa- 
gating within the nucleus are related according to [37|: 

{E - Vf = (Mo + Sf + Pi„2 (5) 

Here V and S denote the real part of the vector and scalar 
potential, respectively. The relation between the momenta 
inside and outside of the nuclear potential are approxi- 
mated by 



M' 



(A/(A/o+^)2+Pi: 



Vf 



(6) 



Refractive effects at the potential boundary were ignored. 
For simplicity no momentum dependence of these poten- 
tials was considered in our schematic simulation. The de- 
fault parameters for the scaler and vector potentials of the 
various baryons at normal nuclear density po are listed in 
Tab. [T] For the proton and the A hyperon (given by 2/3 
of that of the proton) these values give rise to typical total 
potentials of -42 and -28 MeV, respectively. For the an- 
tiproton the summed potential is in the range of the more 
recent experimental results 16|, [iTI, uM, llSl . 



Since the antiproton annihilation and the subsequent 
AA-pair production take place in the nuclear periphery at 
low densities p (see Fig. [T]), the local potentials are ex- 
pected to be reduced. We assumed for simplicity a linear 
density dependence oc p/ po for all vector and scalar po- 
tentials. All numbers for potentials quoted below refer to 
the value at normal nuclear density p^. 



In a last step a finite momentum resolution of 10% was 
applied to mimic possible experimental uncertainties. 

Fig. [3 shows the average transverse momenta of A hy- 
perons (closed symbols) and A hyperons (open symbols) 
as a function of the total A potential V{A) + 5(A) for var- 
ious parameter sets. In all plots the black points result 
from the default parameter set. The transverse momenta 
of the A hyperons drop with decreasing depth of the total 
A potential. For the A hyperons this drop is even more 
pronounced. This surprising behavior can be traced back 
to the imposed momentum conservation and the differ- 
ent sign of the vector potentials for hyperons and antihy- 
perons. Thus within our schematic model an agreement 
between the transverse momenta of hyperons and antihy- 
perons would not necessarily imply that both encounter 
identical potentials. 

Even if all antihyperon and hyperon potentials are set 
to zero one finds different average transverse momenta for 
A's and A's of 225 and 215 MeV/c, respectively. They are 
marked by the dashed and solid arrows in Fig. [2] This dif- 
ference is caused by the anisotropy in the 0cm-distribution 
(Eq. [3|) and the relativistic transformation of the isotropic 
Fermi momentum. We checked that if either the assumed 
Fermi momentum distribution in the target is switched 
off or if an isotropic ^cm-distribution is assumed one ob- 
tains equal average transverse momenta for A and A hy- 
perons. The largest sensitivities to variations of the model 
parameters are observed for the assumed Fermi momen- 
tum (left middle panel) and the assumed anisotropy (left 
lower panel) . For all other parameters the sensitivity of pT 
is significantly weaker. Thus, the simultaneous measure- 
ment of transverse momenta of hyperons and antihyperons 
can be used to adjust the Fermi momentum in the calcu- 
lations. 

Studying only the average transverse momentum dis- 
tributions separately does obviously not allow to extract 
unambiguous information on the potential of antihyper- 
ons. On the other hand, a difference between transverse 
momenta of the concident hadron and antihadron within 
one event reflects directly the different potentials. In or- 
der to study this correlation and to reduce the influence 
of the center of mass angle 9cm we suggest to explore the 
transverse momentum asymmetry ut as a function of the 
longitudinal asymmetry a^. Here, ay and a^ are defined 
for each event in terms of the transverse and longitudinal 
momenta, respectively: 



aT 



Pt(A)-pt(A) 



aL 



Pl{A)~plW .7^ 



PT ( A) + PT ( A) PL (A) + PL (A) 

The Fermi motion of the struck proton inside the target 
may provide a total transverse momentum which will be 
different for each event. However, in case the scalar and 
vector potentials of hyperons and antihyperons are equal, 
the average (ax) is expected to be except for small asym- 
metries caused by the combined effect of the Fermi motion, 
the anisotropic angular distribution and the different ab- 
sorption cross sections for A's and A's. 
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Figure 3: Average transverse momentum asymmetry as a function of 
the longitudinal momentum asymmetry for different parameter pairs 
of tfie scaler and vector A potentials. In each panel calculations with 
3 different Fermi momenta of 180 MeV/c (dashed Unes), 220 MeV/c 
(solid lines), and 260 MeV/c (dotted lines) are overlaid. 



The histograms in fig. [3] show the average transverse 
asymmetry (ax) for various bins in aL- Calculations were 
done for three different scalar potential 5'(A) of -100, -200, 
and -300 MeV, respectively. In each panel, the different 
colored histograms are the results for total potentials of 
V(A) + S(K)^0, -200 and -400 MeV, respectively. Further- 
more, in each plot calculations with three different Fermi 
momenta of 180, 220, and 260 MeV/c are overlaid. For 
all other parameters the default values were used. Varia- 
tions in the first and last bins are partly caused by the low 
number of events in these bins resulting in statistical er- 
rors SaT — 0.01 and 0.02, respectively. While the average 
transverse momenta are very sensitive to the choice of the 
Fermi momentum, the transverse momentum asymmetry 
is not. Using different parametrizations of the Fermi mo- 
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Figure 4: Average transverse momentum asymmetry of AA(solid 
green line), SS (red dashed line) and D+D~ pairs (blue dotted line) 
produced exclusively in 1.66, 2.9 and 6.7 GeV/c p'^'^C interactions, 
respectively. In these calculations the potentials given in the last 
eight columns of Table ^ were used. The black solid, dashed and 
dotted histograms are the result if all scalar and vector potentials of 
the produced AA, SS and D+D~ pairs are set to 0. Asymmetries 
identical to in the first and last bins signal zero counts. 



tion (based e.g. on a local density approximation) gave 
rather similar final results provided the Fermi momentum 
parameter was tuned to similar average transverse hyperon 
and antihyperon momenta. Indeed, for a given A potential, 
ar is mainly determined by the total potential V{A)+S{h.) 
as indicated by the overlap of the histograms of the same 
color. At negative values of a^ the A momenta are rel- 
ative large and consequently the sensitivity of cut to the 
potential is weaker. 

The average transverse asymmetry is non-zero even if 
the total potential V{K) + S'(A)=0 (green histograms). A 
very similar behavior is found even if all scalar and vec- 
tor potentials for A and A are set to 0. Like in the case of 
the average pt, this is caused by the interplay between the 
isotropic Fermi motion and the anisotropic cm-distribution 
and - although less important - by the different absorp- 
tion cross section. Also the assumed scaling of the mo- 
mentum resolution with the absolute momentum causes 
a small positiv correlation between ut and ul- We also 
checked that the results are rather robust despite signif- 
icant changes - by typically ±50% - of all other param- 
eters like absorption cross sections and (anti)proton po- 
tentials. Systematic relative shifts remained usually below 
5a/a = ±0.15. Furthermore, neglecting the momentum 
dependence of the A absorption and assuming constant 
absorption cross sections of lOOmb and 20mb yields asym- 
metries very similar to the momentum dependent cross 
section of Eq. [3] scaled by a factor 1.5 and 0.5, respec- 



tively. 

As already mentioned before, this kinematic method 
of transverse momentum correlations can in principle be 
applied to each hadron-antihadron pair produced exclu- 
sively in pA interactions. As an example, the colored solid, 
dashed and dotted histograms in figure [4] show the average 
transverse momentum asymmetry of A A, SS and D^D~ 
pairs produced in 1.66, 2.9 and 6.7 GeV/c p + ^^C interac- 
tions, respectively. For simplicity, isotropic center-of-niass 
distributions were assumed in case of the SS and D+D~ 
production. For the S and S baryons the same absorption 
cross sections as for the A and A were adopted, whereas for 
D~ and D"*" mesons energy independent absorption cross 
sections of 10 and 90mb, respectively, were taken. The 



scalar and vector potentials were inspired by Refs. |38l. l39l| 
and are listed in Tab. [T] To mimic experimental effects 
a relative resolution for the momentum reconstruction of 
5% was also taken into account. 

For orientation, the black histograms in Fig. 2] show 
the asymmetries if all scalar and vector potentials of the 
outgoing hadrons and antihadrons are set to 0. While for 
SS (dashed) and D+D" pairs (dotted) these histograms 
are symmetric around aL=0, it is not the case for AA pairs 
(solid line). As already mentioned before, this is caused 
by the assumed anisotropic cm. -distribution of the AA 
production. 

Also these calculations confirm the robustness of the 
transverse momentum asymmetry with respect to varia- 
tions of the model parameters: as expected from the sim- 
ilar values for S and V (Tab. [T]), the asymmetries for the 
AA pairs (green histogram in Fig. U]) are indeed close to the 
red histograms shown in Fig. [3] We also would like to note 
that a quite similar result is found for S^S pairs. In a 
purely classical, non-relativistic picture the asymmetry is 
of the order of AC//4 • _Eo, where AC/ is the potential differ- 
ence and Eo the typical kinetic energy of the hadrons. In 
line with this consideration the large laboratory momenta 
of the S hyperons and the D mesons explain the smaller 
asymmetries for the heavier particles. 

To demonstrate the experimental feasibility of the pro- 
posed measurement one may consider as an example the 
bins of 0.25< ar < 0.5 and 0.5< ot < 0.75 where a siz- 
able asymmetry is predicted. Depending on the choice of 
parameters these bins contain approximately 6-10% and 
1-3%, respectively, of the total number of events. The ar- 
distributions have a typical width of 0.3. At the expected 
AA detection rates at PANDA 41] measurement periods 
of a few minutes will be sufficient to reach a relative sta- 
tistical uncertainty of better than 10% for ax within these 
two bins. Given the relative large cross section for pp— *■ 
SS at 2.9 GeV/c of ~l/ib 40] the sensitivity of the trans- 
verse asymmetry (dashed histograms in Fig. |3|) is suffi- 
ciently large to explore the SS pair production at the fu- 
ture FAIR facility. In case of SS pairs a measurement of 
ax with a precision of 10% in the same two bins will re- 
quire typically 2 and 10 hours, respectively. For D-meson 



pairs, however, the large momenta relative to the target 
remnant (> 3 GeV/c) and the low production cross section 
('^lOnb) casts a meaningful measurement of the transverse 
momentum asymmetry in doubt for the case of the poten- 
tial difference of w50-100 MeV given in Tab. [H (cf. the 
dotted black and dotted blue histograms in Fig. |4]). Only 
for significantly deeper potentials than the ones listed in 
Tab.[T]a measurable asymmetry can be expected. But even 
then the estimated measurement periods will significantly 
exceed one month. 

The fact that energy and momentum conservation are 
the main ingredient of the proposed method raises hope 
that similar results might be obtained by more realistic cal- 
culations taking for example the momentum dependence 
of the potentials into account. Since most of the emitted 
hyperon-antihyperon pairs are created in the nuclear pe- 
riphery at subsaturation density, a neutron skin of neutron 
rich target nuclei may help to explore different effective 
potentials. Significant deflections at the potential bound- 
ary which are ignored in the present work may be at least 
partly eliminated by demanding that the target nucleus re- 
mains intact. Furthermore, it may be interesting to study 
questions related to e.g. the formation time [42] by using 
target nuclei of different size. 
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